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Abstract. 

We study numerically the dynamics of a one-electron wave packet in a two- 
dimensional random lattice with long-range correlated diagonal disorder in the presence 
of a uniform electric field. The time-dependent Schrodinger equation is used for this 
purpose. We find that the wave packet displays Bloch-like oscillations associated with 
the appearance of a phase of delocalized states in the strong correlation regime. The 
amplitude of oscillations directly reflects the bandwidth of the phase and allows to 
measure it. The oscillations reveal two main frequencies whose values are determined 
by the structure of the underlying potential in the vicinity of the wave packet maximum. 
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1. Introduction 

Anderson localization theory describes many relevant aspects concerning the nature 
of one-electron states and collective excitations in random media [U El E]. In one- 
dimensional (ID) and two-dimensional (2D) electronic systems with time-reversal 
symmetry, this theory predicts the absence of a disorder-driven metal-insulator 
transition for any degree of uncorrelated disorder. Recently, however, it has been 
reported that the presence of short 01 El El El El CEO] or long-range correlations [HI [121 
[131 [H] in disorder acts towards the appearance of truly delocalized states in ID Anderson 
models. This theoretical predictions was put forward to account for transport properties 
of semiconductor superlattices with intentional short-range correlated disorder [15] and 
microwave transmission spectra of single-mode waveguides with inserted long-range 
correlated scatters fHJ. 

In 2D, the existence of extended states due to correlations in disorder has been also 
proved. Thus, in Ref. [17] the authors considered a two-dimensional striped medium in 
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the (x, y) plane with on-site correlated disorder. The on-site energies were generated by 
a superposition of an uncorrelated term in the x-direction and a long-range correlated 
contribution along ^/-direction. It was predicted that this model displays a disorder- 
driven Kosterlitz-Thouless metal-insulator transition provided strong correlations in 
disorder. More recently, the effect of an isotropic scale-free long-range correlated disorder 
on the one-electron eigenstates of the 2D Anderson model has been studied [18] . To 
introduce long-range correlations in both x and y directions, the site energies were chosen 
to have a power-law spectral density S(k) oc l/k a2D , where k is the magnitude of the 
typical wave- vector characterizing the energy landscape roughness. The metal-insulator 
transition induced by strong correlations (o^d > 2) was monitored by measuring the 
participation number exponent from the long-time behavior of the wave function spatial 
distribution [18] . 

It is well known that a uniform electric field applied to a periodic lattice causes the 
dynamic localization of electron wave packets and gives rise to their oscillatory behavior, 
the so-called Bloch oscillations [TH]. Due to the advances in semiconductor technology, 
it has become possible to monitor the Bloch oscillations in uniform superlattices [20J. 
Remarkably, this phenomenon is not restricted to electronic systems. Recently, the 
authors of Ref. [21] reported the first experimental observation of photonic Bloch 
oscillations in two-dimensional periodic systems. In spite of the fact that the periodicity 
of the potential has been admitted to be the key ingredient for Bloch oscillations to exist, 
it has been demonstrated recently that a ID Anderson model with diagonal long-range 
correlated disorder displays Bloch oscillations in the strong correlation limit |22j. It 
turned out that the period of oscillations agrees well with the one in an ideal Bloch 
band, while the amplitude of oscillations is proportional to the width of the delocalized 
phase, which has been predicted to appear in the strong correlation limit [11]. The 
spectral counterpart of Bloch oscillations — Wannier-Stark quantization of the energy 
spectrum — has been found to be also a remarkable feature of the model [23] . 

In this paper, we focus on the interplay between the derealization effect, arising 
from the long-range correlated disorder, and the dynamic localization, caused by an 
electric field acting on the system. By solving numerically the 2D time-dependent 
Schrodinger equation for the complete Hamiltonian, we compute the behavior of an 
initial Gaussian wave packet in the presence of a uniform electric field. We found clear 
signatures of Bloch-like oscillations [22] of the wave packet between the two mobility 
edges of the phase of delocalized states. The amplitude of the oscillatory motion of the 
centroid allows us to determine the bandwidth of the delocalized phase. 

2. Model and Formalism 

We consider a 2D electron moving in a random long-range correlated potential on a 
regular N x N lattice of unitary spacing and subjected to a uniform electric field F . 
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The corresponding tight-binding Hamiltonian reads [3] 

H = 2_. ( e m + U ■ m) \m){m\ + J 2_, (l m )( n l + \n)(m\) , (1) 

m (mn) 

where \m) is a Wannier state localized at site m = m x e x + m y e y , e m is its energy, 
and U = eF is the energetic bias, e being the electron charge. Here e x and e y are the 
corresponding Cartesian unit vectors. We will assume that the electric field is applied 
along the diagonal of the square lattice. Then U = U(e x + e y )/y2. Transfer integrals 
are restricted to nearest-neighbor sites and are given by J. Hereafter, we fix the energy 
scale by setting J = 1. The long-range correlated sequence of site energies e m are 
generated by making use of a Fourier transform method as follows 

E V2 \^ C(«2d) f2nm x k x {x) \ f2irm y k y {y) \ 

where k 2 = k^ + k 2 , and and 4>^} ky are N 2 /2 independent random phases uniformly 
distributed in the interval [0, 2ir] and C{ a 2D) is a normalization constant, such that 
(e^n) = 1- We also shift the on-site energies in order to have (e m ) = 0. The Wannier 
amplitudes evolve in time according to the time-dependent Schrodinger equation which 
can be written as [3] 

+ U ■ m) -0™ + (VWe* + ^m-e, + VWe s + ^m-ey) , (3) 

with h = 1. Having introduced the model of disorder, we solve numerically Eq. ([3]) to 
study the time evolution of an initially Gaussian wave packet of width a centered at 
site m 

(m — m ) 2 



fp m {t = 0) = A(a) exp 



(4) 



4<x 2 

Once Eq. ([3]) is solved for the initial condition (j4j), we compute the projection of the 
mean position of the wave packet (centroid) along the field direction 

1 

3 



R(t) = - 1 ={e x + e y ) -J2 m \^m(t)\ 2 . (5) 



m 



3. Results and discussion 



In the absence of disorder (e m = 0), all the states are dynamically localized by the 
uniform electric field and the centroid oscillates in time, revealing Bloch oscillations [19|. 
The amplitude and the period of these oscillations are estimated from semiclassical 
arguments (see, e.g., Ref. [24J) as Ly = W/U and tb = 2n/U, respectively, where W 
is the width of the Bloch band in units of the transfer integral J (W = 8 in our case). 
The frequency of the Bloch oscillations is given by u = 2tt/tb and, therefore, is equal 
to the bias magnitude U in the chosen units. 

To study joint effects of the bias and correlated disorder, we performed numerical 
simulations of Eq. ([3]). In all simulations square lattices of size N x N = 250 x 250 were 
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used. The initial wave packet was considered to be located in the center of the lattice, 
m = (N/2)(e x + e y ) and its standard deviation was set to a = 1. 

In Fig. [H (upper panel) we plotted the centroid time behavior of a biased wave 
packet (U = 1) in a weakly correlated random potential (a 2 D — !)• Recall that in the 
absence of bias all the states are localized for such a value of the correlation exponent 
a 2 D [18| - As is seen, switching on the bias does not lead to coherent oscillations of 
the wave packet. A more detailed inspection of these numerical data shows that Bloch 
oscillations in the weak correlation regime (a 2 o < 2) can be observed only for a short 
initial transient. They are strongly damped by disorder and rapidly transformed into an 
incoherent motion of the centroid presented in Fig. [I] (upper panel). To provide a further 
confirmation of this claim, we calculated the Fourier spectrum R(uj) of the centroid R(t). 
Figure [T] (lower panel) shows the result after averaging over 50 realizations of disorder. 
We observe that the Fourier spectrum R{uj) is rather broad, suggesting that R{t) is 
similar to a short-range correlated noise signal with no typical frequency. 
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Figure 1. Upper panel-Centroid dynamics of a biased wave packet (U = 1) in a 
correlated random potential ([2]) with a2D — 1. Lower panel-Fourier transform of the 
centroid, averaged over 50 realization of disorder. 

A rather different time-domain dynamics is found in the limit of strongly correlated 
disorder, a 2 D > 2, when a phase of extended states emerges in the center of the band 
in the unbiased system [18]. We will show below that these states drastically affect the 
dynamics of the wave packet, that now resembles the oscillatory motion of the electron 
in a biased disorder-free lattice. In Fig. [2] and [3] we plotted the centroid time behavior 
of a biased wave packet for two values of the correlation exponent a 2 D = 4 (Fig. [2]) and 
«2D = 5 (Fig. [3]). Two magnitudes of the bias, U = 1 and U = 1.25, were considered for 
each value of a 2 r>. Figures [2] and [3] clearly demonstrate the occurrence of sustainable 
Bloch-like oscillations. Their amplitudes L c are found to be L c w W c /U, where W c is 
independent of the applied bias U. From the data in Fig. [2] we obtain W c ~ 2. This 
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value agrees remarkably well with the width of the band of extended states reported in 
Ref. [IS]. 
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Figure 2. Centroid dynamics of a biased wave packet in a correlated random 
potential ([2|) with ot-iD — 4 calculated for two bias magnitudes U — 1 and U = 1.25 
(left panel). One can see a clear signature of sustainable Bloch-like oscillations. 
The corresponding Fourier spectra of the centroid, obtained by averaging over 50 
realizations of disorder, are shown in the right panel. 
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Figure 3. Same as in Fig. [2] but for a-m — 5. The inset shows the Fourier spectrum 
of the centroid for a single realization, clearly illustrating the fact that the oscillations 
reveal two dominant frequencies around to = U . 

The Fourier spectra R(uj) of the centroid, computed after averaging over 50 
realization of disorder, show a single broad peak at about uj = U . The peak frequency 
nicely agrees with the one in an ideal Bloch band. Looking, however, at the Fourier 
spectrum of a single realization (as shown in the inset in Fig. [3]), one notices that 
actually R(u) has two narrow peaks. The frequency of these peaks fluctuates from one 
realization to the other, which after averaging results in a broad single-peaked spectral 
density. This splitting is not found in ID long-range correlated potentials. 

With the aim to elucidate the anomalies found, we present a simplified model that 
sheds light on the origin of this doublet structure in the Fourier spectrum of the centroid 
oscillations. To this end, we recall that the random site potential e m is given by the 
sum of harmonic terms (T5]). The amplitude of each term decreases upon increasing 
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the harmonic number. For sufficiently high values of «2D, the first term in the series 
will be dominant, while the others are considerably smaller. Consequently, the site 
potential for a given realization represents a harmonic function, perturbed by a colored 
noise (see Ref. [25J). Based on this observation, we keep only the first term in (j2J) and 
neglect all others in our further argumentations. Thus, the "random" site potential now 
is e m = £cos (2%m x /N + <jM') cos (2%m y /N + 0^), where phases 0^ and (j£ v > can be 
arbitrarily chosen. 
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Figure 4. Left panels show the centroid time behavior of a biased wavepacket moving 
in a potential e m — ^cos(2%m x /N + <fi x ) cos(2-Km y / N + <p y ) for the two sets of phases 
(j) x and 4> y indicated in the plots. Right panels show the corresponding centroid Fourier 
spectra. 



The effective local bias is a superposition of the external one and the gradient of 
the local potential, = U — V m e m . At the initial position, m = (N/2)(e x + e y ), 
its components will be given by 

Uf = U+ (2tt(/N) sin (x) cos (y } , (6) 
Uf = U+ (2tt(/N) sin (y } cos (x) . (7) 

For (p^ = cj)( y \ the components of the effective bias are identical. As the frequency 
of the Bloch oscillations is proportional to the magnitude of the effective bias, only a 
single dominant frequency shall be present in this case. One just observes this in Fig. H] 
(left-upper panel). Notice that for the particular case plotted (<^) = 0(f) = 0), the 
local bias at the central position does not acquire any contribution coming from the 
harmonic potential, which results in a null shift of the oscillation frequency. However, 
for the general case of ^ , the local bias components will differ by an amount on 
the order of 1/N. Therefore, Bloch oscillations will have two dominant frequencies. This 
feature is exemplified in Fig. H] (left-lower panel) where we used = and 0^ = 0.5. 
For this case C/J at m is not influenced by the local potential and the corresponding 
oscillation frequency is not shifted. On the other hand, the typical frequency associated 
with oscillations along the y-direction is shifted from the bare frequency u — U = 1. 
The two peak structure depicted in Fig. H] (left-lower panel) and in the inset of Fig. [3] 
(right-upper panel) has, therefore, their origin in the distinct contributions given by the 
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gradient of the local potential to the effective local bias. Averaging over random phases 
(fffl and 4>^ results in a broader Fourier spectrum R(u) around to = U = 1 with a mean 
width on the order of 1/N. 

4. Summary and concluding remarks 

We studied the electron motion on a 2D square lattice with on-site long-range correlated 
disorder in the presence of an external uniform electric field. Long-range correlations 
were introduced by using a 2D discrete Fourier method which generates an appropriate 
disorder distribution with spectral density S(k) oc l/k a ' 2D . Solving numerically the 
Schrodinger equation, the time evolution of an initial Gaussian wave packet was 
investigated, aiming to find coherent Bloch oscillations. Our results suggest that the 
oscillations are strongly damped in the weak correlation limit a^o < 2, when all the 
states are localized because of disorder. The motion of the wave packet on a large time 
scale is chaotic in this case. The sustained oscillations arise for a 2 D > 2, when a phase 
of the extended states emerges at the center of the band. The amplitude of oscillations 
was found to be proportional to the energy difference between the two mobility edges of 
the delocalized phase [18], in good agreement with semiclassical arguments. Thus, we 
arrive at one of the principal conclusions of this work: (i) there exist clear signatures 
of the Bloch oscillations of a biased Gaussian wave packet in the strong correlation 
regime (a2D > 2), originated from the presence of the two mobility edges, and (ii) the 
oscillations exhibit two dominant frequencies because the local bias has a contribution 
from the site energy potential. This is understood using semiclassical arguments. 

The richness of the predicted dynamical behavior can lead to new electro-optical 
devices, based on the coherent motion of confined electrons. We hope that the present 
work will stimulate experimental activities along this direction. 
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